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ABSTRACT 

A review is given of latest results on tests of the flavour independence of strong inter¬ 
actions. Heavy quark mass effects are evident in the data and are now taken into ac¬ 
count at next-to-leading order in QCD perturbation theory. The strong-coupling ratios 
a b s /a^ ds and /a^ ds are found to be consistent with unity. Determinations of the 6-quark 
mass mb(Mz ) are discussed. 
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1. Motivation 


In order for Quantum Chromodynamics (QCD) to be a gauge-invariant renormalisable 
field theory it is required that the strong coupling between quarks (q) and gluons (g), 
a s , be independent of quark flavour. This basic ansatz can be tested directly in e + e“ 
annihilation by measuring the strong coupling in events of the type e + e _ —> qqg for 
specihc quark flavours. Whereas an absolute determination of a s using such a technique 
is limited, primarily by large theoretical uncertainties, to the 5%-level of precision [jlj, 
a much more precise test of the flavour-independence can be made from the ratio of the 
couplings for different quark flavours, in which most experimental errors and theoretical 
uncertainties cancel. Precise measurements of this type have been made using data 
collected at the Z° resonance by the SLC and LEP collaborations, and are reviewed 
here. 

However, the emission of gluon radiation in 66 events is expected to be modified 
relative to that in q{qi (qi =u+d+s ) events due to the large 6-quark mass, and such 
an effect must be taken into account in the test of flavour independence of the strong 
coupling. Recently three groups have completed calculations of 3-jet observables in Z° 
—> bb events complete at next-to-leading order in QCD perturbation theory MSI- 
Comparison of the rates for Z° —> bbg and Z° —y q{qig hence allows measurement of 
the mas^| of the 6-quark at the Z° scale, m^Mz)- Finally, given recent excitement 
in the electroweak sector concerning the values of the quantities R^ and A it seems 
worthwhile to provide a complementary cross-check of the strong dynamics of the 6- 
quark. 

2. Strategy 

The analysis strategy is, in principle, straightforward: one chooses infra-red- and 
collinear-safe observables to measure in 66, cc and q{qi events, and compares the results 
*Use of the modified minimal subtraction renormalisation scheme [^| is implied throughout. 
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with the QCD predictions, taking mass effects into account. However, there are at 
least two possible approaches: 1) input the best knowledge of the 6-mass taken from 
independent measurements, eg. in the T system, and test the flavour-independence of 
the strong coupling; 2) assume the flavour-independence, and determine the 6-mass 
from the data. One or the other of these approaches has been followed in the results 
reported here. 

The basic experimental method is common to both cases. One first selects e + e _ —> 
hadrons events, and applies a flavour-tagging algorithm to select samples of events of 
different primary quark flavour, represented by the indices i and j. One measures the 
3-jet observable(s) X in the different tagged-flavour event samples, and, in order to 
cancel sources of common systematic error, forms the ratios X 1 fX 3 . Before these can be 
compared with the relevant QCD calculations, they must be corrected for: the effects of 
detector acceptance and resolution; the bias of the flavour tag to select preferentially 2- 
jet rather than 3-jet events; the flavour compositions; and hadronisation effects. These 
corrections are very important, and the associated uncertainties must be taken into 
account in the estimation of the systematic errors [|7|, ||. |9]. 

3. Flavour Tagging 

The three new contributions reviewed here have used complementary methods for the 
separation of events of different primary quark flavour, and the systematic uncertain¬ 
ties are therefore at least partly uncorrelated. All methods rely on the large decay 
multiplicity and long lifetime of B hadrons as a discriminator of 66 against cc and qiqi 
events. 

DELPHI has employed a method based on the probability P that all tracks in an 
event originate from the primary e + e _ interaction point (IP) 0; qiqi events tend to 
have a high probability value, and 66 events a low one. Requiring P < 0.005 (P > 0.2) 
tags a 66 (qiqi) event sample with efficiency (e) and purity (n) of 55% and 85% (80% 
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and 80%) respectively. 

OPAL has employed a method based on the number N S i g of ‘significant’ tracks, 
i.e. those whose impact parameter deviates significantly from the IP ||; q{qi events 
tend to have few such tracks, and bb events several. Requiring N sig > 5 ( N sig = 0) tags 
a bb ( q{qi ) event sample with e, II = 23%, 96% (35%, 86%) respectively. 

SLD has employed a method based on the mass, M vtx , and momentum, P vtx , of 
secondary decay vertices reconstructed using the 300M-pixel CCD vertex detector || . 
qiqi events rarely contain reconstructed secondary decay vertices, and these typically 
result from strange particle decays and are of low mass. Conversely, bb events typically 
contain high-mass vertices. Requiring M vtx >1.8 GeV/c 2 tags a bb event sample with 
e, fl = 62%, 90% respectively. Requiring no vertex and N sig = 0 tags a q(qi event 
sample of e, fl = 56%, 91% respectively. 

In addition, both OPAL and SLD have tagged a third sample enriched in primary 
cc events. OPAL has explicitly reconstructed D* candidates with e, II = 2%, 55% re¬ 
spectively. The low efficiency is due mainly to the intrinsically low D* decay branching 
ratios into convenient tagging modes. SLD has utilised the correlation between vertex 
mass and momentum to distinguish cc from bb events by requiring M vtx <1.8 GeV/c 2 
and P vtx > 5 GeV/c. This yields a tagged sample with higher e, II = 19%, 64%, and, 
equally importantly, low bias against 3-jet events. 

4. Test of Flavour Independence 

Having measured the ratios X 1 /X 3 for event flavour samples i and j, and made the 
corrections discussed in Section 2, the flavour-independence of strong interactions can 
be tested via 

V = V(A>; + fl‘«(a!) 2 

Xi AS(X)ai + Bi(X)(ai ) 2 ’ 

where A and B are the perturbatively-calculated leading- (LO) and next-to-leading 
order (NLO) coefficients, respectively, for observable A". In the case of bb (and cc) 
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events these coefficients depend on the quark mass, and, for a given input mass value, 
can be calculated ||. Eq. [I] can then be solved to obtain the ratio a l s /a J s . This has 
been done in new contributions from DELPHI, OPAL and SLD, using Amass values 
of m b (M z ) =2.67 ± 0.27, M b = 5.0 ± 0.5 and m b (M z ) =3.0 ± 0.5 GeV/c 2 , respectively 
(see Section 5). Using the Durham (D) jet-finder with a jet resolution parameter 
value y c = 0.02, DELPHI has measured the 3-jet-rate ratio R™ = R^/R^. Using the 
calculations of Ref. p| they obtained [0] : 

cR 

-f - = 1.007 ± 0.005(exp.) ± 0.007(frag.) ± 0.005(theor.) (2) 


OPAL has measured ratios of the distributions of the event shape obervables ZU, 1 — T, 
Mh , Bt , Bw and C. Using the calculations of Ref. |4j] they obtained 0: 


a; 


a 


uds 


a c 


a 


uds 


= 0.988 ± 0.005(stat.) ± 0.004(syst.) ± 0.011(theor.) 


= 1.002 ± 0.017(stat.) ± 0.025(syst.) ± 0.009(theor.) 


(3) 


(4) 


SLD has measured ratios of the 3-jet rates, i? 3 , for a wide range of y c values, for jets 
defined using the E, E0, P, P0, Durham and Geneva (G) iterative clustering schemes. 
Using the calculations of Ref. 0 they obtained 0: 

,,b 


a" 


a 


uds 


= 1.004 ± 0.018(stat.)lgQ3j > (syst.)io;o 2 9(theor.) 


ac 


= 1.036 ± 0.043(stat.)^o;o45(syst.)lQ;o 2 g(theor. 


oi 


uds 


\ + 0 . 020 / 


(5) 


( 6 ) 


All results are consistent with unity, i.e. no flavour dependence of the strong cou¬ 
pling. The SLD result has larger experimental errors as it is based on a data sample 
comprising approximately 150k Z° events, compared with the roughly 2.8M and 4.4M 
event samples used by DELPHI and OPAL, respectively. The SLD analysis is cur¬ 
rently being updated with the additional 400k events collected between 1996 and 1998 
with the new vertex detector, and the resultant experimental precision is expected to 
be competitive with that achieved by the LEP experiments. Interestingly the three 
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collaborations quote rather different theoretical uncertainties (Eqs. This may 

be due to the facts that DELPHI used only one observable (see next section) and 
implicitly assumed a c s = a^ ds , that different conventions were used for assigning the 
uncertainty, and that different procedures were used by OPAL and SLD for averaging 
over the ensemble of observables. In order to compare results meaningfully and form 
a world average it would be desirable to reach consensus on these issues. The results 
for a b s /a^ ds are illustrated in Fig. |lj. They are also in agreement with older results (not 
shown) incorporating mass effects at LO only [jl]. 


DELPHI 


l l □ l l 

OPAL 

1—1 □ 1 

—I 

SLD 

-1- 

-B-1-1 


0.9 0.95 1 1.05 1.1 

b / uds 

cx s /a s 


Figure 1: Summary of new measurements of a b s /a™ ds determined using quark mass 
effects calculated at NLO. 


5. Heavy Quark Mass Effects 

The effects of a non-zero 6-quark mass are clearly visible in the data. For example, the 
SLD measured ratios R bl are shown in the bottom of Fig. 2(a). For the E, EO, P and 
PO algorithm cases lies significantly above unity (note that all the data points are 
highly correlated with one another), whereas for the D and G algorithms R bl lies just 
below unity. For a 6-mass value of rrib(Mb ) = 3.0 ± 0.5 GeV/c 2 the expected [[J QCD 
values of R% are shown as the short lines below the data points in Fig. 2(a) with the 
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arrows pointing towards lower mass. The calculations are clearly in good agreement 
with the data. Unfolding these data using Eq. || yields the results for a b s /a^ ds shown 
at the bottom of Fig. 2(b) and discussed in the previous section. 


(a) Stat. error only (b) Stat. error only 



Rj/Rf </«f 


Figure 2: R\/R^ ds and a l s /a^ ds measurements for i = c,b. 

One sees, therefore, a 6-mass effect that is jet-algorithm dependent, with R% > 1 
for the JADE family of algorithms and R'f < 1 for the D and G algorithms, and which 
is reproduced beautifully by the NLO calculations. Qualitatively similar effects were 
found for the observables used in the OPAL study || . This can be understood simply 
in terms of two competing physical origins. 1) The non-zero 6-mass tends to cause a 
phase-space suppression of gluon emission relative to the massless quark case. 2) For a 
given kinematic configuration, the large 6-mass tends to enhance the invariant mass of 
a local quark-gluon pair relative to the massless quark case. Since the JADE family of 
jet algorithms is based on a clustering metric that is closely related to invariant mass, 
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for fixed y c the two partons are more likely to be resolved as separate jets when the 
quark is massive, implying R™ > 1. By contrast, the clustering metric used in the 
Durham and Geneva algorithms is less sensitive to this kinematic effect, the phase- 
space suppression dominates, and R% < 1. For increasing values of y c one expects both 
effects 1) and 2) to diminish in importance, and R% —> 1. This is indeed observed, as 
illustrated in Fig. |j from DELPHI, which shows R^ l (y c ) for the D case. 



Figure 3: R% vs. y c compared with QCD calculations using Mb and rrib(Mz )(see text). 


Also shown in Fig. ^jare curves representing the LO and NLO QCD calculations || 
for two choices of the 6-mass based on T data: the pole mass M b (=4.6 GeV/c 2 ) and 
the corresponding running mass m^Mz) (=2.8 GeV/c 2 ). The NLO prediction with 
rrib{Mz) clearly provides the best description of the data, and DELPHI has fitted it to 
the data point at y c = 0.02 by allowing the 6-mass mb{M z ) to vary. They obtained 

mb(M z ) = 2.67 ± 0.25(exp.) ± 0.34(frag.) ± 0.27(theor.) GeV/c 2 


(7) 











This result is illustrated in Fig. [|, where it is plotted at the energy scale /i = Mz and 
compared with the value [}7j m^Mx/ 2 ) = 4.16 ± 0.14 GeV/c 2 . The comparison yields: 

rrib(Mz) — mb(Mx/ 2 ) = —1.49 ± 0.52 GeV/c 2 (8) 

which is consistent with the expected QCD running, shown by the band in Fig. (|. 



Figure 4: The running of the 6-mass with energy scale. 


DELPHI has updated its study recently 0 with an improved flavour tag, and 
has also employed the new Cambridge jet-finding algorithm 0- The corresponding 
dependence of on y c is shown in Fig. || In this case, although the NLO calculation 
with mb(Mz) describes the data best, the NLO and running mass effects do not seem 
as pronounced as in the Durham case, and both the LO calculation with rrib(Mz) and 
the NLO calculation with M h are also consistent with the data. 

The running 6-mass has also been studied 0| using the SLD data shown in Fig. 2. 
The calculations described in Ref. [fl have been performed, for the six jet algorithms 
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and y c values used, for mb(M z ) values in the range 2.0 < rrib(Mz) < 4.0 GeV/c 2 . The 
results are illustrated in Fig. |^, where the curves are cubic polynomial interpolations 
between the calculations at discrete mass values. It can be seen that the rrib(Mz)- 
dependence has a different form for each algorithm, with a positive slope for the E, 
E0, P and P0 cases, and a negative slope for the D and G cases. By comparing these 
curves with the SLD data, shown as horizontal bands, for each algorithm one can read 
off the value of rrib(Mz) that is preferred, represented by the vertical lines. 



Figure 5: vs. y c compared with QCD calculations using M h and mb(M z )( see text). 


It is worth reemphasising that the SLD results shown are highly correlated among 
the six jet algorithms since the same data sample was used and the observables are 
intrinsically correlated. It is therefore interesting that no single 6-mass value fits the 
data corresponding to all six algorithms. The r.m.s. scatter of the central values 
is Arrib(Mz) = ±0.49 GeV/c 2 , which is comparable with the DELPHI total error on 
m,b(Mz ) using just the Durham algorithm. The origin of this scatter is unclear, but it is 
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certainly plausible that it results from higher-order QCD effects that are, by definition, 
not included in the NLO calculation, and which can a priori be of different sign and 
magnitude for different obsrvables. In addition, the considerable dependence || on 
jet algorithm of the size of the hadronisation uncertainties on R b J translates [12] into 
uncertainties on rrib(M z ) as large as GeV/c 2 . 

A weighted average over the results from the six jet algorithms yields a preliminary 
result: 


m b (M z ) = 3.23lo;72( s tat.)l?;28(syst.)15;of(theor.) ± 0.49(r.m.s.) GeV/c 2 (9) 

This is consistent with Eq. |7j. and the r.m.s. has been conservatively included as 
an additional uncertainty. (For this reason the r.m.s. over the six a b s /a^ ds results 
discussed in the previous section was included as an uncertainty on the SLD average 
value, Eq. |). Finally, it should be noted that effects due to the c-quark mass are 
expected to be much smaller than those due to the Amass, i.e. at the level of 1% or 
less on i?g/ R ^ ds . Such effects are much smaller than the current experimental errors, 
and are difficult to observe; see Fig. 2 and Refs. |8|, |9|]. 


6. Summary and Conclusions 

Recent results on tests of the flavour independence of strong interactions and heavy 
quark mass effects from DELPHI, OPAL and SLD have been reviewed. Three theoret¬ 
ical groups have recently completed calculations of 3-jet observables, including quark 
mass effects, complete at next-to-leading order in perturbative QCD. These calculations 
have been used enthusiastically by the experimentalists in the flavour-independence 
studies, as well as to determine the Aquark mass at the Z G scale. 

The values of the strong coupling ratios, a b /a™ ds and a. c s /a.™ ds , extracted using the 
different NLO calculations, are in good agreement with one another, with the previous 
results incorporating mass effects only at LO, and with unity. There is hence no 
evidence for any anomalous flavour-dependent effects. Though the measurements of 
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Figure 6: SLD R\/R^ ds measurements compared with the mfe(M^)-dependence of the 
NLO calculation (see text). 
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a b s /a™ ds are approaching the per cent level of precision, there are fewer measurements 
of a c s /a ™ ds , and these are limited in precision currently to the 5% level; more and better 
measurements would be welcomed. It would also be desirable to achieve a consensus on 
the treatment of experimental and theoretical systematic uncertainties by the different 
collaborations, so that the data could be combined sensibly to obtain meaningful world 
average values and errors. 

The DELPHI and SLD data have been used to determine m,b(Mz ) arid consistent 
values are obtained. A strong jet-algorithm dependence of mb(M z ) has been observed, 
as well as large hadronisation uncertainties for some jet algorithms. DELPHI is cur¬ 
rently investigating use of the new Cambridge jet algorithm, and the recent large SLD 
data sample collected with the new vertex detector will be included in an analysis 
optimised for the study of rrib(M z ). 
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